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We develop a general variational approach to study the statistical properties of the tail states of 
a wide class of non-Hermitian operators. The utility of the method, which is a refinement of the 
instanton approach introduced by Zittartz and Langer, is illustrated in detail by reference to the 
problem of a quantum particle propagating in an imaginary scalar potential. 
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Over recent years considerable interest has been shown 
in the spectral properties of random Fokker-Planck op- 
erators, and their application to the dynamics of vari- 
ous classical systems ||l|,|[. Of these, perhaps the best 
known example is the "Passive Scalar" problem which 
concerns the diffusion of a classical particle subjected 
to a random velocity field. Here, in contrast to quan- 
tum mechanical evolution, the random classical dynam- 
ics is typically specified by a linear non-Hermitian oper- 
ator. Non-Hermitian operators also appear in a number 
of problems in statistical physics. For example, the sta- 
tistical mechanics of a repulsive polymer chain can be 
described in terms of the classical diffusion of a parti- 
cle subject to a random imaginary scalar potential [^,||. 
Similarly, the statistical mechanics of flux lines in a type 
n superconductor pinned by a background of impurities 
can be described as the quantum evolution of a particle in 
a disordered environment subject to an imaginary vector 
potential M. While these connections have been known 
for a long time, the ramifications of non-Hermiticity on 
the nature of the dynamics is only now being fully ex- 
plored. Beginning with early work on random matrix en- 
sembles [pl0l, a number of attempts have been made to 
analyze spectral properties of non-Hermitian operators, 
and to apply the results to the description of classical sys- 
tems |5|,|8|-[l9| . For example, it has been found that the 
generic localization properties of non-Hermitian systems 
are drastically different from those of their Hermitian 
counterparts, a fact that can be attributed to an implicit 
chiral symmetry of non-Hermitian operators |^,^,g^ . 

Previous studies of non-Hermitian operators have been 
largely based on perturbative schemes, such as the self- 
consistent Born approximation in the diagrammatic anal- 
ysis (e.g. Ref. [|l3| ), or the mean field approximation 
in the field-theoretic approach (e.g. Ref. Il^l)- How- 
ever, there are indications 12^-03] that an important role 
can be played by those parts of the spectrum which are 
populated by exponentially rare, localized states. These 
"Lifshitz tail" states, first introduced in the context of 
semiconductor physics, can not be treated perturbatively. 
Since the original works of Lifshitz pj| , a number of so- 
phisticated mathematical methods to deal with tail states 
have appeared pSpq] , including the instanton technique 
in statistical field theory pTl. The aim of this letter is 



to propose a new non-perturbative scheme to investigate 
properties of the tail states of non-Hermitian operators. 
Although our approach is quite general, for clarity, we 
choose to explain the main features of the method by ap- 
plying it to possibly the simplest model system. Specif- 
ically, we study the Hamiltonian of a quantum particle 
subject to a random imaginary scalar potential 



H = -A + iV{r) 



(1) 



where the potential V{r) is drawn from a Gaussian 6- 
correlated impurity distribution with zero average, and 
correlator given by {V{r)V{r')) — 'yS{r — r'). Previous 
studies have shown that, when averaged over realizations 
of V{r), the Feynman propagator of the Hamiltonian 
above dl^) can be identified with the partition function 
of a self-repelling polymer chain with a contact interac- 
tion HM]. Moreover, this model can be used to describe 
NMR in inhomogeneous materials, where the interplay of 
diffusion and local variations in the spin precession rate 
results in the dynamics of magnetization being specified 
by the operator (nl) ES] . Indeed, in the case of the lat- 
ter, the tails of the operator (H) can be shown to provide 
a significant contribution to the relaxation of the NMR 
signal. 
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FIG. 1. Support of the states in the complex plane {d — 1). 
The density of states outside the shaded region is zero if cal- 
culated in the self-consistent Born approximation. 

The Hamiltonian (1^) is non-Hermitian, and its eigen- 
values Efe occupy some area in the complex plane. In the 
self-consistent Born approximation, the density of com- 
plex eigenvalues, defined as p = (X]fe ^(^ ~ ^k)6{y — yt)), 
is given by [R^ 



p(x,y) 



I (47r7)- 
10, 



l2/|<A(x), 
\y\ > A(a^), 



(2) 



where x and y are the real and imaginary parts of the 
complex energy e — x + iy, A(a;) = 2Trju{x), and u^x) 
is the density of states of the clean (7 — 0) system. The 
validity of this "mean field" result is restricted to values 
of X such that A{x) <?C x, which corresponds to large x 
in dimensions lower than four. 

According to Eq. ^, there are no states outside the 
region \y\ < A (a;). This can be understood if one recalls 
that the mean field approximation is essentially a way 
to determine the self-consistent contribution of typical 
fluctuations of the random potential. On the other hand, 
numerical simulations Q suggest the existence of rare 
states with \y\ ^ A(a;), which are generated by atypically 
strong fluctuations of V{r). Such fluctuations occur with 
an exponentially small probability exp[— M^/7], where 



W^^ I V^{r)dr. 



(3) 



For any given configuration V{r) of the random poten- 
tial, we will describe W as the energy associated with 
that configuration. This should not be confused with the 
"energy" e of a quantum state, which by definition is the 
eigenvalue of H corresponding to that state. Since strong 
fluctuations of the disorder potential are suppressed by 
an exponential factor, states with I2/I S> A(x) are domi- 
nated by those configurations of V{r) that have the high- 
est statistical weight exp[— 14^/7] or, equivalently, the low- 
est energy W. Thus, we come to the idea of the opti- 
mal fluctuation method, whose rigorous formulation can 
be given as follows: Amongst all the configurations of 
V{r) such that e is an eigenvalue of the Hamiltonian P), 
we choose the one that minimizes the energy functional 
VF[F(r)]. The density of states is then given with expo- 
nential accuracy by the statistical weight of the minimal 
configuration, p{e) ^ exp[— Wmm(e)/7]. This approach 
is similar to that employed in the seminal work by Zit- 
tartz and Langer |29] in the treatment of the Hermitian 
model. There, a saddle-point technique was applied to es- 
timate the functional integral over V{r). Although based 
on the same physical ideas, this technique leads to incon- 
sistent results when formally applied to a non-Hermitian 
problem. Instead, we introduce below a more general 
variational formulation tailored to the consideration of 
non-Hermitian operators. 

The arguments outlined above are not restricted to the 
particular form of the Hamiltonian p), and can be ap- 
plied whenever one has to deal with tail states (i.e. local- 
ized states created by strong fluctuations of the disorder 
potential). Formally, the limitations of this scheme are 
set by the inequality W„iin{x,y) 2> 7, which defines a 
certain area in the complex plane of eigenvalues. Below 
we will show that the precise form of the optimal fluctua- 
tion potential can, in general, be determined by solving a 



system of coupled non-linear equations. For the particu- 
lar case of the Hamiltonian (nh, we will obtain an explicit 
form of the solution in the limit y <^ x. 

With this preparation, we now turn to the derivation 
of the variational approach: Our aim is to minimize the 
functional T1^[F(r)] subject to the constraint det(i?— e) = 
0. Account for the latter can be made by introducing two 
Lagrange multipliers, /ii and ^2, which in turn leads to 
the functional 

F[V{v)] = W - niRe det{H - e) - /i2 Imdet(ff - e). 

(4) 

If e is an eigenvalue of H with left and right eigenfunc- 
tions tp^{Y) and ■(/'^(r), a spectral decomposition of the 
Hamiltonian (p]) obtains the identity 



det(iJ-e)=*V^^*(r)^''(r) \{{ek-e), (5) 



5V{r) 



k,<Lk=^e 



where efc denote the remaining eigenvalues of H (i.e. 
those different from e). Equating the functional deriva- 
tive 5F/5V{v) to zero, and applying Eq. (|^) we obtain 

y(r) = AiRei/'^*(r)V'^(r) + A2lmV'^*(r)V'^(r), (6) 

where Ai and A2 denote redefined Lagrangian multipliers. 
When combined with the eigenvalue equations 

HV''^(r) = eV^(r), i?V/(r) = e* '0^(r), (7) 

and with the binormality condition 

drV'^*(r)V'^(r) = l, (8) 



Eq. (|^) defines a closed system of non-linear equations for 
-0^, ^■^, Ai, A2, which has to be solved in order to extract 
the optimal configuration ^(r) of the disorder potential. 
Note that Eq. (||) depends explicitly on the nature of 
the impurity distribution and on the structure of the op- 
erator at hand, whereas Eqs. (0) and (H) are universal. To 
treat other random non-Hermitian operators, one simply 
has to modify Eq. (o) accordingly. For example, in the 
case of the Passive Scalar problem, where the disorder 
has the form of a random velocity field, the analogue of 
Eq. (ra) relates the optimal configuration of the random 
velocity field to spatial derivatives of the eigenfunctions 

For the Hamiltonian (nh, further progress can be made 
by exploiting symmetry properties. From the relation 
W = H* it follows that one can constrain the complex 
wavefunctions to obey the relation -(/^^(r) = 0'^*(r) — 
4>{^)- In doing so, leaving aside the lengthy but straight- 
forward analysis |30|, one can show that the following 
simplifications obtain: Firstly, Eq. (H) can be effectively 
disregarded; secondly, by performing a complex rotation 



of tp{v), one finds that arbitrary non-zero values can be 
assigned to the Lagrange multipliers Ai and A2. For con- 
venience, we will set Ai = 1, A2 = 0. With these simpli- 
fications, focusing initially on the one-dimensional prob- 
lem, one finds 



V = Reip^, 



(9) 



where the wavefunction ip is obtained self-consistently by 
solving the Schrodinger equation in the potential V , 



ip" + iVtp = eip. 



(10) 



Eqs. ^ and (|l(]|), which represent the non-Hermitian 
analogue of the non-linear Schrodinger equation of Zit- 
tartz and Langer, are the main result of the paper. To 
complete the program, one should find the localized so- 
lution ip (which is assumed to be unique) of Eqs. (|^) and 
( pl3| ) for each value of e, and calculate the energy W{e) of 
the corresponding configuration V — Re ip'^ of the poten- 
tial. Note that, in the derivation of Eqs. (0) and (ttO), no 
approximation has been made. However, the applicabil- 
ity of the optimal fluctuation method itself is restricted 
by the condition Wmi„(e) » 7. 

As mentioned above, an analogous calculation can be 
carried out for other model systems including the Passive 
Scalar operator 



Hps = -A 



-v(r) -v. 



(11) 



Indeed, in the case of a Gaussian distributed incompress- 
ible (V • V = 0) flow, the tail states of Tips are found to 
be governed by equations which have the same form as 
Eqs. (I) and (0) 0. 

In the analysis of Eqs. (^) and (pi)|), it is conve- 
nient to interpret the one-dimensional spatial coordi- 
nate r as a time t, and the complex wavefunction tp{t) 
as the position of a fictitious classical particle in the 
two-dimensional plane. With this interpretation, one 
can recast Ec^s. <M) and (nO) in the Lagrangian form, 
idC/dilj[ 2)' = dC/diJi,2, where 



C = lj)[%j}'2 - xV'lV'2 



|(.? 



^1) 



i(^? 



r^Y, (12) 



and tjj ~ tjji + i^j2 has been separated into real and imag- 
inary parts. The Lagrangian (12) has at least one invari- 
ant of the motion — the classical energy. Although we 
can not rule out the possibility of this system being in- 
tegrable, we have been unable to find a second invariant 
of the motion. Therefore, integrability remains an open 
question, which deserves a separate investigation. 

As follows from Eqs. (pi) and (Eoh, the wavefunction of 
a tail state is defined by the components of the energy, 
X and y, which determine the position of the state in 
the complex plane of eigenvalues. In order to exploit the 
optimal fluctuation method to its fullest potential, one 
should analyze these equations over the entire complex 



plane. However, at present, we do not have a method of 
finding the exact analytical solution of Eqs. (^) and ( [lO| ) 
for arbitrary values of x and y. Therefore, we will focus 
on a specific domain, in which an approximate solution 
can be found. 

As we will see later, the density of Lifshitz tails de- 
cays exponentially as a function of the distance from 
the boundary of the "mean field spectrum", \y\ = /S.{x). 
Therefore, the most physically relevant states lie close to 
the boundary. For such states, the condition \y\ -^ x 
is satisfied (recall that A(a;) <C x) allowing the following 
separation of scales: Since the wavefunction of a tail state 
oscillates with a high frequency ^Jx, while the amplitude 
of these oscillations varies at a much lower rate \y\/ ^/x 
(the ratio of the two time-scales being \y\lx <^ 1), the 
solution of Eqs. (^ and (|l^) can be parameterized in the 
form of a wavepacket. 



^(t) 



<^H 






^t^t ^ ^_ 



Uvij " 



■i\fxt 



(13) 



Substituting the parameterization (O) into Eqs. (S) 
and (|l^), one obtains an equation for the envelope. Fol- 
lowing a standard procedure ^^, we neglect the sec- 
ond derivatives of (/j± , and perform the averaging over 
the intermediate scales. As a result, one obtains a sys- 
tem of parameter-free first-order differential equations for 
V±(t): 






- ip+ + 5tp\ip- + ipt- 



f- - ^P+V 



^^ 



(14) 



In general, the solution of equations of this kind, which 
can easily be computed numerically, yields a universal 
dimensionless constant in the final result for the density 
of states. Surprisingly, for the problem at hand, Eqs. ( |l4| ) 
can be represented in the form of Hamilton's equations, 
wherein ip± play the role of the canonical variables, (p'j. — 
ztdH/difizp, with the Hamiltonian 

H - -^+^- + ^ifilipl + ^{ifiX + ipi). (15) 

Of the infinite set of solutions, defined by values of 7i, we 
are interested in the localized solution, which corresponds 
to Ti = 0. Integrating Eqs. ( p^ along the line ?i = 0, we 
obtain 



'/'±(t) = 



2e 



±T 



Vi 



At 



+ e- 



-At 



10 



(16) 



Now, with the help of Eq. (O) , we can restore the wave- 
function tp{t). It has a rather unusual shape shown in 
Fig. 2(a). The combination of an oscillatory wavefunc- 
tion modulated by a localized envelope is characteristic 
of the tail states of non-Hermitian operators. The op- 
timal configuration of V(t), which is expressed through 



tpit) according to Eq. (^, also has the form of locahzed 
oscillations (see Fig. 2(b)), and its energy can be straight- 
forwardly determined: 

VK™„(e) = a\y\V^, " = 4f ^^^(^ + 2\/6) « 3.744. 

(17) 

Recalling that, in one dimension, ^{x) = l/(27rY^), we 
arrive at the expression for the density of tail states. 



p{e) r^ exp 



\y\ 



A(x) 



A(a;) < \y\ < x. 



(18) 



We remark that, in higher dimensions, the spherical 
symmetry of the optimal potential reduces the problem 
to being effectively one-dimensional. In this case, one 
obtains the general result pQ] , 




p\v\ \y\ 



X A(a;) 



A(a;) < |y| < a;, 

|y|>2;, 

(19) 



where the values of the numerical factors a and /3 depend 
on the dimensionality, and A(a;) = 2-k^v{x) ^ ^x^'^~'^^l'^ . 
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FIG. 2. Solution of the non-linear Scrodinger equation (g) 
and ( [L0| ) for a; = 8.0, j/ = 1.0: (a) shows the square modulus 
li/^P of the complex wavefunction. The optimal fluctuation of 
the disorder potential V = Re-i/)^ is shown in (b). 

For completeness, one should add that in addition to 
the tail states depicted in Fig. 1 and studied here, the 
random imaginary scalar potential exhibits another class 
of exponentially rare states which inhabit the region of 
very small a; ^- ||I^ . The origin of the latter seems to be 
associated with large areas in real space which are almost 
free of disorder. Such states are inaccessible within the 
present framework. 

To summarize, in this paper we have developed a gen- 
eral variational approach to study spectral properties of 
non-Hermitian operators. Applied to the problem of a 
particle propagating in a random imaginary scalar po- 
tential, one finds that properties of the strongly localized 
(tail) states are governed by a system of coupled non- 
linear dynamical equations. These equations show that 



tail states associated with the non-Hermitian Hamilto- 
nian exhibit oscillations on scales much shorter than the 
localization length. Employing a procedure of scale sep- 
aration, both the density of tail states as well as the 
corresponding wavefunctions have been determined. 

We would like to thank V. Ruban, K. Samokhin, 
L Smolyarenko and A. Moroz for helpful discussions. 
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